Electronic structure of a single vortex in d-wave superconductor revisited 
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The present work deals with the study of d-wave superconductor in presence of a single vortex 
placed at the centre of the 2D lattice using the t — t' — J model within the renormalized mean field 
theory. It is found that the in absence of the vortex the ground state has a d-wave configuration. 
In presence of the vortex, the superconducting order parameter, above the critical doping, drops 
to a low non zero value within a few lattice points from the vortex (that is within the the vortex 
core) and beyond it converges to the constant value in absence of vortex. We observe that above 
the critical doping things are consistent with the experimental results while there is an anomalous 
rise in the superconducting order parameter at very low doping within the vortex core. This we feel 
is because of antiferromagnetic ordering taking place within the vortex core at low doping. 
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Introduction: In a type-II superconductor, magnetic 
field enters into the system in terms of vortices. The or- 
der parameter owing to broken symmetry in a supercon- 
ductor aquires a topological phase of 2tt for an winding. 
The structure of the vortices in superconductors is tra- 
ditionally an interesting subject of research in condensed 
matter physics. In a conventional BCS superconductor 
with s-wave pairing symmetry [l| , the formation of quasi- 
particle bound states inside the core of a vortex that had 
been predicted long back, has been verified in scanning 
tunneling microscopy (STM) measurement 

But a similar convincing picture in an unconventional 
superconductor like cuprate superconductor with d x i_ y i- 
wave pairing is still elusive. This is because the pair 
potential in the latter kind of superconductor is non- 
local and it has nodes. The conventional weak cou- 
pling BCS theory for d-wave vortices predicts zero en- 
ergy peak in the local density of states at the vortex 
core for the Bogoliubov quasiparticles from the solution 
of the Bogoliubov-de Gennes (BdG) equations 6]. How- 
ever, the scanning tunneling microscopy measurements 
on the cuprate superconductors seem to reveal a bound 
state inside the vortex core at fairly high energy 0, Q ■ 
This fact is counter intuitive, since quasiparticle exci- 
tation is possible even at zero energy in a pure d-wave 
superconductor 0,J1] . There have been several proposals 



like eL 



state, a s tagg ered flux state, and an 
antiferromagnetic vortex core [lOjto describe the absence 
of ZEP inside the core. Subsequently, neutron scattering 
experiments near optimal doping ll| , ttSR and NMR ex- 
periments in underdoped compounds [12j seem to suggest 
the existence of the enhancement of the AF correlation 
inside the core [t| Therefore the absence of zero 

bias tunneling conductance in the STM measurements 
may be due to the AF gap developed inside the core. In 
the present work we study the d-wave superconductiv- 
ity using the t — if — J model using renormalized mean 
field theory (RMFT). In this scheme we use Gutzwillcr 
projected wavefunction which renormalizes the parame- 



ters of the Hamiltonian. The resulting Hamiltonian is 
decoupled using the Unrestricted Hartree-Fock scheme 
and then the corresponding Bogoliubov-DeGennes equa- 
tions are solved in a self consistent fashion. The central 
idea is to capture the strong electronic correlation effi- 
ciently. We study the simplest case to begin with where 
decoupling is done only in favour of the superconduct- 
ing order parameter. Depending on the results of this 
study we intend to include the decoupling in favour of 
antiferromagnetic order in future in the required regime 
of doping. 

The Model: We consider a Hamiltonian 

i,($i,er i,<52,er i,<7 
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with no doubly-occupied sites. Here c\ a represents cre- 
ation of an electron with spin a =\ or \, at site i, 
Si = ±x, ±y, 82 = ±i ± y. (Here the inter-site distance 
is chosen to be unity.) The hopping energies t and t' are 
for the nearest neighbor (NN) and next nearest neighbour 
(NNN) respectively, and fi is the chemical potential which 
fixes the average density of electrons n = 1 — x with hole 
doping x. The spin operator Sf = \ J2 a /3 Aa^a^P witn 
o~ a being the Pauli matrices, and the density opcartor 
7ii = c,Lc,;t + cJiCj.. The exchange energy J term 



in Eq. ([T]) may be decomposed into NN bond pair- 
ing amplitude | [(c^c i+(5li )o - (cjj.Cj+^o] which de- 
scribes d x 2_ y 2 -superconductor with a mean field state 
$o)- The suppression of double occupancy is achieved 
by the Gutzwiller projection of |$o) into the state 
1$) = Hj(l — rii^ni±)\&o} ■ The expectation values in 
this new state is related with that of ordinary mean field 
state as follows: ($14^ |<I>) = g t ($ \4a c 3<r\®o) and 
<*|Si • S 3 \$) = ff,(# |St: £4$o>, where g t = 2x/(l + x) 
and g s = 4/(1 + x) 2 HI]. Therefore the Guzwiller 
renormalized mean field Hamiltonian for a superconduc- 
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tor becomes 

ft = c L c i+si* + t' E c I CT ^+5 2 ^ -aE c ^ Cict 
+E [ A *.i+*i { c lt4+8u - 4A + s lt ) + h - c ] > ( 2 ) 
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where i — gti^ i — gtt\ ft is the renormalized chemi- 
cal potential which fixes x for Ji and the effective pair 
amplitude 
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with J = (3p s + l)/4. The Hamiltonian (@J can be diago- 
nalized using standard Bogoliubov transformatin leading 
to the usual BdG equations for both kind of Bogoliubov 
quasiparticles: 
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where E n is the energy for the n-th eigen value, and 
w" and v™ are the corresponding amplitudes at the i-th 
site for creation of a spin-f electron and destruction of 
a spin- 1 hole respectively. The NN pairing potential can 
be obtained through the Gutzwiller meanfield equation 
as 



[E n + fl)u n (¥^G. 1: Configuration of a vortex placed at the center of 
a 40 x 40 lattice. Black dot represents the position of the 
centre of the vortex. Length and angle of an arrow represent 
magnitude and phase of the pair amplitude respectively 
at the centre, which is denoted by the tail of the arrow, of 
two neighbouring lattice site. 
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For a pure d x i_ y i superconductor, Ai t i±$ = — A^j±g = 
Ao, a constant. If we introduce a vortex, pair amplitude 
aquires a phase in the center of mass coordinate (CMC) 
of the pair: A iri ±% = A e l6, ' i±s and A iti ±y = -Aoe'^ 1 ** 
where fty is the angle of the CMC of i and j with respect 
to the center of the vortex. This implies phase winding of 
2tt around the vortex. Choosing realistic parameters as 
t' = t/4 and J = i/3, we then solve the BDG equations 
for eigenvalues E n and eigenfunctions {u l n ,v l n ) for differ- 
ent values of mu in the presence of a vortex placed at the 
center of an N x N (N even) lattice. Even N is chosen 
for placing the vortex in the center of an atomic lattice 
as it is the most symmetric and energetically favorable 
position. We begin with the above form of A^^ for a 
vortex and solve the BDG equations in the open bounday 
condition together with self-consistency for the average 
hole concentration x = -As — ^ J2 n K 
erage NN pair amplitude A = 2jv( ^_ 1) 6i 
have iV = 40, 50, and 60 in our finite size computation. 
However, there is no boundary effect on the phase as has 
been shown in Fig. [T] for N=40. There is no boundary 
twist on the phase. The bulk magnitude of pair ampli- 
tude is independent of N for the chosen numbers. This 
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2: The behavior of A(r) in the unit of its value at the 
as a function of the distance r away from the center 
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FIG. 
bulk 

of the vortex for different values of x. The values of x are 
0.05, 0.08, 0.13, 0.31, and 0.49 from top to bottom at the 
left edge of the curves. The inset shows A(r) when y — 0.5 
for N= 40 (circle), 50 (square), and 60 (triangle) at the hole 
concentration x « 0.05 



is demonstrated in the inset of Fig. [2j We therefore be- 
lieve that we have achieved the thermodynamic limit to 
predict the electronic structure of a single vortex in high 
T c superconductors. 

For the largest system with N = 60 that we have con- 
sidered, the magnitude and phase of the pair amplitude 
between the electrons in the bonds around a single vortex 
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FIG. 3: Structure of a single vortex in a 60 x 60 lattice with 
hole concentrations x = 0.05, 0.13, 0.31, and 0.49 (from top 
to bottom panels). Black dot represents the center of the 
vortex. Blue dots are the position of the Cu atoms in the 
lattice. Length, direction, and tail of the arrows represent 
the magnitude of bond-pair amplitude, phase of bond-pair 
amplitude, and the center of the bonds respectively. 
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FIG. 4: The bulk values of pair-amplitude as a function of 
x (upper-curve). Its renormalization with the factor g t = 
2x/(l + x) is shown in the lower-curve. 



is shown in Fig. [3] for the hole concentrations 0.05, 0.13, 
0.31, and 0.49. The phase winding around the vortex in 
each case is 2ir. The magnitude of the pair amplitude 
A(r) is constant away from the core of the vortex, i.e., 
A(r — > oo) = A . Its dependence on x is shown in Fig.0] 
However, the behavior of A(r) inside the core have sev- 
eral features depending on x. (i) At the highly under- 
doped regime, unlike any other superconductor, A(r) is 
greater than Ao at the closest bond from the center of the 
vortex, it is lesser than Ao in the immediate next bond 
and the next bond onwards it is almost close to Ao. (ii) 
Near quantum critical point, x = 1/8, A(r) is almost in- 
dependent of r. (iii) Near optimal doping, the behavior 
of A(r) inside the core is exactly opposite to the case 
of highly underdoped regime, (iv) At the highly over- 
doped regime, its behaviour near the center of the vortex 
is the same as conventional vortex, i.e., its value gradu- 
ally decreases towards the center of the vortex. Since the 
calculation is performed in a lattice, the value of A(r) 
will not be the exactly same for same values of r but at 
different bonds. We therefore perform best fit for A(r) 
and have demonstrated its behaviour as a function of 
r for different values of x in Fig. [21 A(r) is scaled by 
Ao for different values of x. The following remarkable 
features are found. (i)The structure of the core at highly 
overdoped regime is conventional, i.e., the pair amplitude 
approaches zero towards the center of the vortex and then 
smoothly attains its asymptotic value, (ii) In highly un- 
derdoped regime, the pair amplitude has larger value at 
deep inside the core compared to the same in the bulk 
and it oscillates about the bulk value before attaining the 
asymptotic value, (ii) The regime around optimal doping 
shows that has smaller value (but does not converge to 
zero as r —> 0) at deep inside the core compared to the 
same in the bulk and it oscillates about the bulk value be- 
fore attaining the asymptotic value, (iv) Near quantum 
critical point in the underdoped regime, pair amplitude 
remains almost constant everywhere. From the curves 
of Fig. 2, we estimate the coherence length to be about 
3 times atomic lattice constant. This is consistent with 
the experiments. Figure HI depicts the behaviour of the 
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FIG. 5: Quasiparticle density of states as a function of quasi- 
particle energy for x = 0.05, 0.13, 0.31, and 0.49 respectively 
from bottom to top seen at E = 0. 

absolute value of the (i-wave order parameter away from 
the vortex core as a function of the doping x with (lower 
curve) and without (upper curve) the renormalizing fac- 
tor being multiplied with it. We observe that the dome 
like feature in the d-wave order parameter is obtained 
with the renormalizing factor g t = 2x/(l + x) being mul- 
tiplied as shown in the lower curve. This is again at par 
with the experimental results. So, as far as the bulk be- 
haviour of the d-wave parameter is concerned the present 
scheme of decoupling the Hamiltonian only in favour of 
the superconducting order parameter is sufficient. The 
antifcrromagnetic ordering, we feel is important for the 
underdoped regime to explain the physics inside the vor- 
tex core. The quasiparticle local density of states (DOS) 
is given by: 

= [KI 2(S ( e " E n) + KMe + E n )] (7) 

n 

Figure [5] show the DOS at one of the four nearest 
atomic site from the center of the vortex for N=60, and 
x = 0.05, 0.13, 0.31 and 0.49. For highly overdoped 
regime, the DOS at the core has large ZEP as was found 
by the conventional study. Near optimally doped regime, 
the DOS is large but does not show any peak near zero 
energy. The underdoped regime shows that the DOS is 
vanishingly small at zero energy but shows a peak at 
E Kt 0.03t. This peak is consistent with the large tunnel- 
ing current observed at fairly high energy. 

Conclusion 

The present work reveals that above the critical dop- 
ing, decoupling the Hamiltonian only in favour of the 
superconducting order parameter is giving results con- 
sistent with the experiments for both the bulk and the 



vortex core. Below critical doping the behaviour of the 
d-wave order parameter in the bulk is consistent with the 
experimental findings but inside the vortex core it shows 
anomalous behaviour as it rises sharply on approaching 
the vortex centre. We believe that in the underdoped 
regime the antifcrromagnetic ordering is important to 
get the correct physics inside the vortex core. Hence 
decoupling the Hamiltonian both in favour of the anti- 
ferromagnetic and superconducting order parameter and 
then solving it self consistently would give results at par 
with the experiments in the vortex core. 
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